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o : Abstract 

^ [ In recent papers we have introduced a method for the study of hmit cycles of the Lienard 

1^ ■ system : x = y — F{x) , y = —x , where F{x) is an odd polynomial. The method gives a 

O I sequence of polynomials Rn{x), whose roots are related to the number and location of the 

■ limit cycles, and a sequence of algebraic approximations to the bifurcation set of the system. 

In this paper, we present a variant of the method that gives very important qualitative and 
^ ■ quantitative improvements. 
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In a previous paper we have introduced a method for studying the number and 
location of limit cycles of the Lienard system : 

dx . dy 

-dt=y-^^''^ ' rf^ = -^ ' 

where F{x) is an odd polynomial of arbitrary degree. The method is as follows : we 
consider a function hn{x,y) given by : 

hn{x, y) = y'' + gn~i,n{x)y'^~^ + 5'„_2,n(a;)y"~^ + ... + gi,n{x)y + go,n{x) , (2) 

where gj^n{x), with j = 0,l,...,n — 1, are functions of x only and n is an even integer. 
It is always possible to choose the functions gj,n{x) such that : 

d • dh dh 

■^hn{x, y) = K{x, y) = {y- F{x))-^ - x-^ (3) 

is a function of the variable x only (see also 0). Hence we have : 

hn{x,y) = Rn{x) (4) 

The functions gj^n{x) and Rn{x) determined in this way are polynomials. We have 
shown in [Q] and that the polynomials hn{x,y) and Rn{x) give a lot of information 
about the number and location of the limit cycles of @). In particular, we have 
established in [jl| the following conjecture : 

Conjecture : Let L be the number of limit cycles of (|l|). Let r„ be the number of 
positive roots of i?„ (with n even) of odd multiplicity. Then we have : 

• L <rn 

• ii m> n then — r„ = 2p with p E N 

Moreover, we have also shown in |Jl|] and that the polynomials hn{x,y) and Rn{x) 
allow us to construct algebraic approximations to each limit cycle and to the bifurca- 
tion curves of ([^). For the bifurcation set, these algebraic approximations are exact 
lower bounds and seem to converge in a monotonous way to the exact bifurcation set 
of the system. The fundamental aspect of this method is that it is not perturbative in 
nature. It is not necessary to have a small or a large parameter in order to apply it. 

In the present paper, we want to improve the results presented in and [Q]. Let 
us consider, as a first example, the van der Pol system : 

X = y — e(x^/3 — x) 

y = -X (5) 
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In this case we have F{x) = e(x^/3 — x) and the polynomials e) have only one 

positive root of odd multiplicity for each even value of n and for arbitrary values of 
the parameter e (we have indicated here the explicit dependence of i?„ in e by writing 
Rn{x,e)). We call amplitude of the limit cycle, the maximum value of the coordinate 
X on the limit cycle and we will refer to it by Xmax- For the van der Pol equation, this 
amphtude is a function of e and we will write : Xmax{^) 

In table (|l|), we give the roots of the polynomial Rn{x,e) for n between 2 and 20 
and for e = 3. These sequence of roots seems to converge in a monotonous way to 
the amplitude of the limit cycle Xmax{^) — 2.023, which is obtained by a numerical 
integration of the system. 

As it is explained in the integral of each polynomial Rn{x,e) along the limit 
cycle must be zero for all even values of n : 

/ Rn{x{t),e)dt = , (6) 

where T is the period of the limit cycle. 

Let us now describe for this example the improved method that represents the new 
contribution of this paper. We employ an idea utilized in the averaging method @] : 
we replace x{t) by a ■ cos(t) in (^), where a is an unknown constant, and we replace 
the period T by 2tt. After integration, we obtain a polynomial in a which we denote 
Rn{a,e) : 

Rn{a,e)= / Rn{a ■ cos{t), e)dt (7) 
Jo 

Surprisingly enough, the polynomials Rn{a,e) have the same qualitative properties as 
the polynomials Rn{x, e). Each of them has only one positive root of odd multiplicity 
for arbitrary values of e. The values of these roots, for n between 2 and 20 are given 
in table (Q) for the case e = 3. We can verify that each one of these roots represents a 
lower bound for x^axi^ = 3). This sequence of roots seems to converge to x^axi^ = 3) 
much more rapidly than the sequence of roots of the polynomials Rn{x,e) and they 
represent excellent approximations to the value of the amplitude of the limit cycle. 

This behavior of the roots of the polynomials Rn{a,e) is merely an experimental 
fact. At present, we have no rigorous arguments to explain these results. We have 
observed this behavior of the roots of the polynomials Rn{o) for other Lienard systems 
of type d^) and the conjecture established in |l[ (and given also above) about the roots 
of the polynomials Rn{x) seems to be valid also for the roots of the " averaged " 
polynomials Rn{ci). 
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For a given value of n, we can obtain the approximation of the amphtude Xmax{^) 
as a function of e by considering the curve given by the equation Rn{x,e) = 0. How- 
ever, a better approximation is found by considering _R„(a, e) = instead. In fig.(0) 
(respectively (@)), we give the curve -R„(x,e) = (respectively i?„(a,e)=0) for several 
values of n and the numerical curve Xmax{^) obtained from a numerical integration of 
the system. As we can see from these figures, the improvement obtained with the new 
method is very important and has two different aspects : 

• a qualitative aspect : the curves -R„(a, e) = are nearer to the numerical curve 
than the curve -R„(x, e) = 0. 

• a qualitative aspect : the asymptotic behavior (when e — oo or e ^ 0) of the 
curves Rn{x, e) = is the correct one (even for small n). 

The amplitude Xmaxi^) of the limit cycle of the van der Pol equation tend to the value 
2 when e ^ oo or e ^ : 

lim Xmax{e) = \imxmax{e) = 2 

This asymptotic behavior is correctly given by the curves Rn{a,e) = for all even 
values of n. 

Despite the fact that the curves Rn{x,e) = do not have the correct asymptotic 
behavior, each one represent an exact lower bound to the function x^axi^) and is closer 
to it than its predecessor. Moreover, for a given value of e, if we take n sufficiently 
large , the root of -R„(x, e) = will be as near as we want to Xmaxi^)- 

For other recent results about the limit cycle of the van der Pol equation see [Q and [^]. 

Let us consider a second example : 

X = y — e {x^ — \/ax^ + x) 

y = -X (8) 

where e and a are arbitrary positive parameters. This system has been carefully studied 
by Rychkov and can have at most two limit cycles. Since there are two parameters, 
the bifurcation set is given by a curve in the parameter plane (e, a). 

In we have shown that the method presented in allows us to obtain a sequence 
of exact algebraic lower bounds to the bifurcation set of systems like (H). Here, we will 
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show that by using the polynomials Rn{a), instead of the polynomials Rn{x), we can 
considerably improve the results presented in [§. 

In the first quadrant of the plane (e, a), there exists a bifurcation curve B{e, a) = 0. 
On this curve, the system undergoes a saddle-node bifurcation (see for a description 
of this type of bifurcation). Obviously, this function B{e, a) = is not known and no 
analytical method for obtaining it for arbitrary values of the parameters exists. 

We will obtain algebraic approximations to the curve B{e, a) = from the polyno- 
mials Rn{x, e, a) and Rn{ci, e, a). We will call Bn{e, a) = the algebraic approximations 
obtained from the polynomials Rn{x,e,a) and 5„(e,a) = the curves obtained from 
the polynomials -R„(a, e, a). As explained in [Q, the function -B„(e, a) is obtained from 
the conditions : 

dR 

Rn{x, e,a) = , —r^{x, e, a) = (9) 
dx 

In the same way, the function -B„(e, a) is obtained from the conditions : 

dR 

Rn{a,e,a) = 0, —-^{a,e,a) = (10) 
dx 

The algebraic equations @ (respectively (p!0D ) determine the double root of R 
(respectively Rn{ci, e, a)) and give a relation between e and a, which we write Bn{€, a) = 
(respectively Bn{e,a) = 0). The curves -B„(e,a) = are shown in fig.(|D for several 
values of n. The curve B{e, a) = 0, calculated from numerical integration of the sys- 
tem, is also given. In figure (^, we show the curves i?„(e, a) = and B{e,a). Again, 
the improvement obtained with the polynomials i?„(a,e, a) is very important. The 
curves i?„(e, a) = represent better approximations to the curve B{e, a) = than the 
curves Bn{e, a) = do. 

It can be proved by perturbation methods that the asymptotic behavior of the func- 
tion B{e, a) when e ^ oo is given by B{e, a) ~ a — 5. The curves Bn{e, a) = have not 
this asymptotic behavior when e oo. On the contrary, the curves -B„(e, a) = that 
we have studied {n between 2 and 20) have a correct asymptotic behavior (see fig.(^). 
In this way, the curves Bn{e, a) = have the right global shape when compared to the 
numerical bifurcation curve. Both families of curves i?„(e,a) = and i?„(e,a) = 
give lower bounds to the unknown exact bifurcation curve B{e, a) = 0. For the curves 
Bn{e, a) = 0, this result has been established in 0. But for the curves Bn{e, a) = it 
is an experimental fact that cannot be proved in a simple way. 

Let us point out that, despite the fact that the curves -Bn(e, a) = have not the correct 
asymptotic behavior , for a given value of a, the value of e obtained from the equation 
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i?„(e, a) = can be as near as we want to the exact value of the bifurcation curve 
provided that we take n sufficiently large. 

In summary, the curves S„(e, a) — represent a sequence of algebraic approxima- 
tions to the bifurcation curve B{e, a) — 0. These approximations are very good, even 
for small values of n. They are better than the approximations given by the exact lower 
bounds Bn{€, a) — 0. The improvement obtained from the polynomials Rn{(i, e, ck) is 
very surprising because it seems that it is a general fact, valid for arbitrary odd poly- 
nomials F{x). The mathematical justification of this method (more specifically the 
passage from the polynomials Rn{x,e,a) to the "averaged" polynomials i?„(a,e, a)) 
represent an interesting open problem. 

In the meantime, the method presented in this paper gives a very effective way 
for obtaining information about the number of limit cycles, their amplitudes and their 
bifurcations for the Lienard systems. 
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n 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 


num. 


root 


1.732 


1.819 


1.863 


1.890 


1.909 


1.923 


1.934 


1.943 


1.950 


1.955 


2.023 



Table 1: Values of the roots of Rn{x, e) for system @ with e = 3 



n 
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10 
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14 


16 


18 


20 


num. 


root 


2 


2 


2.003 


2.006 


2.008 


2.010 


2.011 


2.012 


2.013 


2.014 


2.023 



Table 2: Values of the roots of -R„(a, e) for system (|^) with e = 3 
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1 . 95 
1 . 9 

1 . 85 
1 . 8 

1 .75 
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Figure 1: Plots of Rn{x, e) = for n = 2 to n = 16 for system (^). The point line is 
the Xmaxle) calculated numericaly. 
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Figure 2: Plots of Rn{cL, e) = for n = 2 to n = 16 for system (^. The point line is 
the Xmax(e) calculated numericaly. Note that, here, the horizontal axe is a = 2. 
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Figure 3: Plots of Bn{(y,e) = for n = 2,6 and 10 for system (|^). The point line is 
B{a,e) = calculated numericaly. 
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Figure 4: Plots of Bn{(y,e) = for n 
B{a,e) = calculated numericaly. 



2,6 and 10 for system (|^). The point line is 
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